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ABSTRACT 

This paper investigates the steady-state responses of a ro- 
tor system supported by auxiliary bearings in which there 
is a clearance between the rotor and the inner race of the 
bearing. A simulation model based upon the rotor of a pro- 
duction jet engine is developed and its steady-state behavior 
is explored over a wide range of operating conditions for var- 
ious parametric configurations. Specifically, the influence of 
rotor imbalance, support stiffness and damping is studied. 
It is found that imbalance may change the rotor responses 
dramatically in terms of frequency contents at certain oper- 
ating speeds. Subharmonic responses of 2nd order through 
10th order are all observed except the 9th order. Chaotic 
phenomenon is also observed. Jump phenomena (or double- 
valued responses) of both hard-spring type and soft-spring 
type are shown to occur at low operating speeds for sys- 
tems with low auxiliary bearing damping or large clearance 
even with relatively small imbalance. The effect of friction 
between the shaft and the inner race of the bearing is also 
discussed. 

NOMENCLATURE 

Cb = auxiliary bearing support damping, lb.s 2 /in. 

Cb%p = auxiliary bearing torsional damping, Ib.in.s 
F n = normal force, lb 
F t = friction force, lb 

Fx = external force vector acting on the rotor in X direction 
Fy = external force vector acting on the rotor in Y direction 
I a = rotor inertia matrix 


Jb = moment of inertia of auxiliary bearing, Ib.in.s 2 
Kb = auxiliary bearing support stiffness, lb/in. 

Kc = contact stiffness, lb/in. 

Mb = auxiliary bearing mass, lb.s 2 /in. 

M k = mass of kth rotor element, lb.s 2 /in. 

N = total number of modes considered 
NB 1 = node number at auxiliary bearing #1 
N B2 = node number at auxiliary bearing #2 
Qx = rotor modal coordinate vector in X direction 
Qy = rotor modal coordinate vector in Y direction 
Rb = radius of auxiliary bearing bore, in. 

Rm = radius of auxiliary bearing pitch, in. 

Rb = radius of rotor journal, in. 

Xfl = rotor physical coordinate vector in X direction 
Yh = rotor physical coordinate vector in Y direction 
e as rotor imbalance eccentricity, in. 
g = gravitational acceleration, in./s 2 
t = time, s 

A = deformation at the contact point, in. 

r = * T i a * 

= rotor free-free modal rotation matrix 
fi = rotor operating speed, rad/s 
$ = rotor free-free modal displacement matrix 
<5 ss R b — Rrj auxiliary bearing clearance, in. 
fi = dynamic friction coefficient 
= rolling friction coefficient 

\ 1>b = angular displacement of auxiliary bearing mner-race 
£ = modal damping coefficient 
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INTRODUCTION 

One of the most innovative developments in the turbo- 
machinery field involves the use of active magnetic bearings 
(AMB) for rotor support. This technology provides the po- 
tential for significant improvements in the dynamic behavior 
of rotor systems, allowing for loading, eccentricity, shaft po- 
sition and vibration to be continuously monitored and con- 
trolled. In order to protect the soft iron cores of the magnetic 
bearings and to provide rotor support in the event of failure 
of the bearing or during an overload situation, backup (or 
auxiliary) bearings, with a clearance between the rotor and 
the inner race of the bearing, are usually included in the ro- 
tor design. This clearance introduces a nonlinear dynamical 
feature which may significantly impact the behavior of the 
rotor. 

Magnetic bearing systems appear to provide particularly 
great promise for use in aerospace applications. There are 
active programs at many of the major jet engine manufac- 
turers to develop engines supported by magnetic bearings. 
Safety is a major concerns in any aeronautical design. To- 
ward this end, it is desirable to design the rotor system to 
take maximum advantage of the backup bearings and use 
them as true auxiliary bearings to provide support during 
critical situations in a safe and consistent manner. An im- 
portant concern in this regard is the dynamic behavior of the 
rotor when it comes into contact with the auxiliary bearing. 
If safe and effective operation of the engine is to be ensured 
during these periods, it is essential that designers have a very 
good understanding of the steady-state dynamics of rotor 
systems with clearance effects. 

There are a number of studies in the literature concerned 
with the dynamics of rotors with clearance effects. Ya- 
mamoto (1954) conducted a systematic study of rotor re- 
sponses involving bearing clearance effects. Black (1968) 
studied the rotor/stator interaction with a clearance. He 
concluded that rotor/stator interactions may occur in a vari- 
ety of forms and circumstances, including jump phenomena. 
Ehrich (1966) reported the first identification of a second 
order subharmonic vibration phenomenon in a rotor system 
associated with bearing clearance (1966). Bently (1974) pub- 
lished experimental observations of second and third order 
subharmonic vibration in a rotor system. Later, Muszynska 
(1984) cited the occurrence of second, third, and fourth order 
subharmonic responses in a rotor rubbing case and Ehrich 
(1988 and 1991) observed eighth and ninth order subhar- 
monic vibration as well as chaotic vibration in a high speed 
turbomachine. Childs (1979 and 1982) published two papers 
to explain the mechanism for the second and third order 
subharmonic responses noted above. He stated, with great 
insight, that “ motion due to nonsymmetric clearance effects 
is a fractional-frequency phenomenon.” 

While those studies have greatly enhanced the understand- 
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FIG. 1 DIAGRAM OF THE FEM ROTOR MODEL 

ing of clearance effects on rotor dynamics, a more detailed 
understanding of the dynamical behavior of such systems is 
needed. The perspective of much of this earlier work is that 
the clearance exists as a result of manufacturing error or mis- 
fitting. That is, it is due to an abnormal situation. However, 
in a rotor system fitted with magnetic bearings and auxiliary 
bearings, the clearance becomes a design parameter rather 
than an irregularity. From this point of view, it is impor- 
tant to develop a detailed quantitative understanding of the 
dynamic responses that are to be expected. Such knowledge 
will provide guidelines for the selection of auxiliary bearing 
parameters. 

It seems that there have been little work to date that is 
specifically concerned with auxiliary bearings in magnetic 
bearing supported rotor systems. Two papers that are di- 
rectly related to research on auxiliary bearings were both fo- 
cused on transient responses. Gelin et al. (1990) studied the 
transient dynamic behavior of rotors on auxiliary bearings 
during the coast down. Ishii and Kirk (1991) investigated 
the transient responses of a flexible rotor during the rotor 
drop after the magnetic bearings become inactive. In both 
papers, idealized rotor models are used and it is assumed 
that once the magnetic bearings fail, the torque is cut off 
and consequently the rotor speed approaches zero. 

In this paper, simulation results are presented for a 
complex rotor system supported by auxiliary bearings with 
clearance at each end of the rotor. This work is specifically 
concerned with systems in which the clearances are quite 
small (on the order of a few mils), which is appropriate 
for jet engine applications in which the backup bearing is 
acting to provide rotor support on a consistent basis. The 
influence of rotor imbalance, support stiffness and support 
damping are investigated using direct numerical integration 
of the governing equations of motion and the harmonic 
balance method. Some insights are obtained with regard 
to the frequency and amplitude behavior of the steady-state 
vibration of such a system. 



2 



coordinate transformation: 



FIG. 2 AUXILIARY BEARING MODEL 

SIMULATION MODEL 

The rotor is modelled using free-free normal mode shapes 
and natural frequencies obtained through finite element 
analysis. The model data is representative for the rotor 
of a jet engine. Fig. 1 shows a schematic diagram of the 
FEM rotor model. Parametric information about the model 
is listed in Table 1. The torsional motion of the shaft is not 
considered in this paper. Using state space representation 
and modal coordinates, the equations of motion for the rotor 
are expressed as 


N 

Xfuc = ^2 Qxi, 

i = 1 
N 

YRk = ^2 Qn* 

»'ss 1 


(Jfc = NBl,NB2) 


The equations of motion for the auxiliary bearings are 
derived using the model shown in Fig. 2 


M B kX Bk *F CakXQk + KakXak 

= F nk cos a k - F tk sin a k + M Bk g, (2. a) 
MskYBk + CakYBk -F X Bk Y Bk 

= F nk sin a k + F tk cos a t , (2.6) 

JBk^Bk + C Bx j, k rp Bk 

= F\ k R Bk fi^F nkRmky (2.c) 


where 


at = tan “ l 


Y Bk - Y Bk 
X Bk — X B k 


{k = NBl y NB2) 


At this point, the rotor and the back-up bearings appear 
to be uncoupled. However, the force vectors and F y 
on the right hand sides of equations (1) are partially due to 
rotor/auxiliary bearing interaction. In fact, we have 


Qx+2Cw n Q* -f QTQy + w nQx 
+ 2f)C^nQy = ^ r Fx, 
Qy-f2Cu> n Qy — 

-2n(u n Qx = ^Fy, 


where 

F* = {Fxi,Fxi, • ■ • , •f’xm} -1 , 

Fy = {FyuFyi fy m } -1 , 

Qx = 

Qy =$- 1 Y fi> 

with 

X/l = {Xfli, Xri, . . . , .YiJrn} -1 , 

(m = total number of nodes) 


(i.«) 

( 1 . 6 ) 


Fxk = - F nk cos a k + F tk sin a k 
+ Mk g + Mk eU 2 cos (fi<), 

Fy k = - Fnk sinatk - F tk cosa k 
4- Mk e Q 2 sin (fll). 

The rotor/bearing interaction is represented with the normal 
force F n k 


where 


Fnk = 


r K C 6k, 

10. 


A* <0, 
A fc >0, 


(3 a) 


A* = (Xjik — XBk) cos oik + (YRk — YBk) sin ak — 6k 


and the Coulomb friction force F tk . As long as there exists 
slip at the contact point, the friction force obeys 


F t k = » Fnk. (3.6) 


The physical displacements of the rotor at the two aux- 
iliary bearing locations can be obtained using the following 
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However, when there is no slip at the contact point, the 
friction forces are solved from equations (1) and (2) using 



the kinematic constraint that the circumferential velocities 
of the rotor and the inner-race of the back-up bearing at 
the contact point equal to each other. At the same time, if 
this solved friction force exceeds the maximum static friction 
force (= fi s Fnk)* equation (3.b) applys again. 

DISCUSSION OF RESULTS 

The rotor is modeled with 34 stations (as shown in Fig. 1) 
and the first four modes (two rigid body and two flexible 
modes) are included in the simulation model. The two 
auxiliary bearings are located at nodes 3 and 33, respectively. 
This arrangement is taken to represent one of the most 
technically feasible configurations in that it greatly simplifies 
bearing maintenance. It is assumed that the two auxiliary 
bearings are identical in terms of stiffness, damping and 
friction characteristics. Some nominal system parameters 
used for the simulation study are /fc=‘2.855e-b6, £=0.03, 
Rm jfc — 1 • 1 Rski Vs— 0-5 and ^=0.002. To avoid excessive 
cluttering of plots, all the results that are presented in this 
paper correspond to node 3, the location of bearing #1. 

Since the total system which includes two bearings and 
associated friction forces as well as the inner-race motions 
is rather complicated and requires considerable amount of 
computer time for the solutions to converge, the friction 
effect is examined first to see if the model can be further 
simplified. It turns out that the steady-state results obtained 
with and without friction are virtually identical. Even the 
differences in transient responses are quite small, as can be 
seen in Fig. 3(a) and 3(b). The only remarkable effect is 
on the transient responses of the inner-race as shown in 
Fig. 3(c). This observation is confirmed by numerous runs 
using different system parameters and rotor speeds. The 
lack of significance of friction may be attributed to several 
factors. First, the inertia of the inner-race is quite small in 
comparison to the rotor mass, so vibration of the bearing 
has little influence on the rotor vibration. Second, ball 
bearings exhibit quite negligible torsional resistance under 
normal conditions. As a result, the terms that are related 
to the friction forces and the rotational motion of the inner 
race are not included in the simulations that are discussed in 
the following paragraphs. 

The steady-state response characteristics of the system 
are obtained through numerical integration of the simplified 
version of governing equations (1) and (2). Near-zero initial 
conditions are used, simulating situations where the AMBs 
are functioning properly prior to a system failure. Multiple 
solutions with other initial conditions are not sought at 
present. 

It is well known from linear analyses that imbalance 
greatly affects the steady state vibration amplitudes of a 
rotor system. However, it is observed from the current 
work that imbalance may also influence frequency content of 


the rotor responses quite dramatically at certain operating 
speeds. A typical case with such an imbalance effect is shown 
in Figs. 4 and 5, where orbits and corresponding frequency 
spectra of the rotor for different values of imbalance at the 
speed of Q = 1000 are plotted. For this particular case, 
there exist eight ranges of imbalance values that result in 
eight different types of rotor responses. 

For e < 0.0009, the rotor rotates near the bottom of 
auxiliary bearings and the responses are predominantly syn- 
chronous (Figs. 4(a)-4(b) and 5(a)-5(b)). As imbalance in- 
creases, the 2Q superharmonic component approaches the 
magnitude order of the synchronous component (Figs. 4(b) 
and 5(b)). However, the responses are of small amplitude. 
For 0.0010 < e < 0.0013, the responses are dominated by 
Q/2 subharmonic components (Figs. 4(c) and 5(c)). In other 
words, the amplitude of Q/2 component is greater than that 
of the synchronous. For 0.0014 < e < 0.0016, the Q/2 
subharmonics disappear and the Q/3 subharmonics become 
dominant (Figs. 4(d) and 5(d)). So far, the overall amplitude 
of the responses are not large, the rotor just bounces near the 
bottom of the auxiliary bearings. For 0.0017 < e < 0.0027, 
the orbits become chaotic-looking (Fig. 4(e)) and the spec- 
trum contains a lot of noise (Fig. 5(e)). In this range of 
imbalance, the rotor changes from bouncing near the bot- 
tom to bouncing around the full clearance of the bearing as 
imbalance increases. In the middle of this transition range, 
true chaos is observed. The Poincare map shown in Fig. 6(a) 
and the frequency spectrum shown in Fig. 6(b) demonstrate 
that the response has all the characteristics of a chaotic 
phenomenon. It should be noted that even though the or- 
bits are chaotic looking, the amplitudes are not the largest 
among all the cases for this particular parametric configu- 
ration. For 0.0028 < e < 0.0034, the orbits orbits are no 
longer chaotic-looking (Fig. 4(f)). The spectrum shows they 
are Q/5 subharmonic responses (Fig. 5(f)). Notice the am- 
plitudes are the largest for this parametric configuration. For 
e = 0.0035, the amplitude suddenly becomes smaller even 
though the imbalance has become larger (Fig. 4(g)). And 
the rotor bounces near the bottom of the auxiliary bearings 
again. The frequency spectrum shows it is Q/8 subharmonic 
response (Fig. 5(g)). For 0.0036 < e < 0.0042, the orbits 
become chaotic-looking again (Fig. 4(h)). But the frequency 
spectra are very similar to the subharmonic cases (Fig. 5(h)), 
only with some discrete noise. Finally, for e > 0.0043, the re- 
sponses become predominantly synchronous again (Figs. 4(i)- 
4(j) and 5(i)-5(j)). But this time as imbalance increases, the 
2Q superharmonic component become smaller and smaller 
(Figs. 4(j) and 5(j)). 

Examining all the orbits in terms of amplitudes as imbal- 
ance increases, we can see the characteristics of a jump-type 
phenomenon (Cunningham, 1958). The jump-down takes 
place around 0.0034 < e < 0035 where the rotor jumps 
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from full-clearance bouncing to near-bottom bouncing. Fur- 
ther investigation is needed to better understand this type 
of change. 

Imbalance responses at some other operating speeds and 
for other parametric configurations exhibit similar changes 
as imbalance varies, though the imbalance ranges and cor- 
responding response types may not be as well defined as 
in the above cases. In fact, subharmonic responses from 
Q/2 through Q/10 are all observed except Q/9 as shown 
in Figs. 7. Surprisingly, those subharmonics are not di- 
rectly related to the system’s natural frequencies as were the 
cases with other researchers’ findings (such as Ehrich, 1988). 
Moreover, several types of subharmonic responses may occur 
at a single operating speed. It should be noted that Chen et 
al. (1993) also reported occurrence of three stable subhar- 
monic responses at a single rotor speed in a SFD supported 
rotor system but did not provide any explanation for their 
findings. In their case, even the imbalance did not vary. Ap- 
parently, further research is needed to find a mechanism to 
explain this multi-subharmonic vibration phenomenon. On 
the other hand, it should be pointed out that these subhar- 
monic responses are not typical cases. While some of them 
are observed to exist within a certain range of parameters, 
the majority of them occur only for some specific parametric 
configurations. 

Due to space limitations, results for other parametric 
configurations are not systematically plotted. A general 
summation of the observations is presented instead. A 
common feature among all the responses is that for very 
small imbalance, the responses are always synchronous. 
The imbalance range that result in synchronous dominated 
responses depend on several system parameters. For small 
back-up bearing stiffness (such as Kq = 0.213e-|-6) and 
normal damping (Cb = 157.0), the responses are almost 
always synchronous. Only Q/2 subharmonic are observed 
at a few operating speeds with a very narrow range of 
imbalance. It should be noted that even though a lower 
Kb may leads to a better system response, it may also 
fail to protect the magnetic bearings due to the fact that 
it could result in a larger rotor orbit-center offset. However, 
the dramatic response changes discussed above may occur 
again if the damping becomes small (such as Cq = 57.0) 
even though the stiffness still remains small. On the other 
hand, increasing the damping Cb alone may not be able 
to eliminate those dramatic changes. It is observed that 
those changes can still occur for Cb being as large as 
700.0. Reducing the size of clearance 8 may not eliminate 
the response changes at certain speeds. But it can narrow 
the operating speed range where those changes occur. For 
example, response changes are eliminated for Q > 1500 
when 6 is reduced from 0.002 to 0.001 with all other 
parameters remaining the same, but response changes still 


occur for Q < 1400. 

It is obvious that rotor responses involving nonsymmetric 
bearing clearance effect are very complex problems and 
numerical integration alone is not a sufficient tool to obtain 
a global picture of the system responses. The harmonic 
balance method is then used for the investigation of global 
system behavior. However, it is only attempted for situations 
with very small imbalance values. In addition, only the lQ 
harmonic is considered. The complex frequency contents 
associated with medium and large imbalance values makes 
it a formidable task to apply the harmonic balance method 
for other cases. Nevertheless, some useful information can be 
drawn from these results. After all, an adequately balanced 
rotor system should have very small imbalance under normal 
conditions. 

Fig. 8(a) shows that nonsymmetric clearance effect is 
equivalent to asymmetric stiffness effect with regards to 
critical speeds. The clearance actually splits the first critical 
speed into two pseudo— critical speeds. In the X direction, 
the gravity force tends to keep the rotor in contact with 
the bearing at low operating speed. Thus, the apparent 
stiffness is almost the same as Kb and th e pseudo-critical 
speed is nearly the same as the critical speed for the linear 
case (8 ~ 0). But in the Y direction, the clearance 
results in a lower apparent stiffness and, consequently, an 
additional lower-value pseudo-critical speed. It is seen that 
several higher order additional pseudo-critical speeds are 
created in the operating speed range in addition to the 
1st additional pseudo-critical speed. It is noted that the 
response in the X direction also departs from the linear 
case at high operating speed. This is because the imbalance 
force becomes dominant at high rotor speed which in turn 
makes the gravity force less significant and the clearance 
effect more important. Fig. 8(b) shows that changing the 
auxiliary bearing stiffness has little effect on the pseudo- 
critical speeds of the system. However, for a larger value 
of imbalance, a higher Kb does leads to a greater tendency 
of double-valued responses. In each direction, for either a 
stiffness increase or an imbalance increase, the 1st pseudo- 
critical peak tend to become a hard-spring type jump and 
the 2nd one tend to develop into a soft-spring type jump, 
with the tendency decreasing as the pseudo-critical’s order 
increases. 

Fig. 9(a) shows double— valued responses in the Y direction 
for four different values of clearance. It is seen that a larger 
clearance results in wider rotor speed range of double-valued 
responses. It is also observed that as clearance increases, 
the apparent stiffness decreases and the first pseudo— critical 
speed shifts to a lower value. Fig. 9(b) shows the double- 
valued responses in the X direction. Even though the jumps 
themselves are smaller in magnitude, they are more obvious 
in trend. Notice how little the change is for the first pseudo- 
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critical speed in the X direction. 

Fig. 10(a) shows the influence of auxiliary bearing damping 
on the double-valued responses in the Y direction. It is 
observed that the damping has to be quite large to eliminate 
the double-valued responses associated with the first pseudo- 
critical speed. Fig. 10(b) shows the influence of Cb on the 
double-valued responses in the X direction. In both figures, 
it should be noted that as Cb decreases, the second pseudo- 
critical speed peak will develop into a soft-spring type jump 
and the third pseudo-critical speed peak will evolve into a 
hard-spring type jump. 

The system behaviors for higher operating speed range are 
not shown in Figs. (9) and (10) so that the jump phenomena 
can be more clearly illustrated. It is also because that the 
system’s responses at high operating speed range with the 
same parameters are more or less regular, in other words, 
mainly amplitude changes. 

CONCLUSIONS 

As a summary of the results discussed above, the following 
conclusions can be drawn: 

1. Imbalance may change the rotor responses dramatically in 
terms of frequency contents at certain operating speeds, 
especially under conditions of large clearance, high bearing 
stiffness and low bearing damping. With imbalance 
changing, as many as eight different types of responses 
may occur for a particular parametric configuration at a 
single operating speed. 

2. Subharmonic responses of second order through tenth 
order are all observed except for the ninth order case. 
However, the majority of them are not typical cases, 
and were observed only for quite particular parametric 
configurations. 

3. Chaotic phenomenon is observed to occur occasionally. 
However, the amplitudes associated with such motion are 
not among the largest. 

4. Nonsymmetric clearance effects influence the critical speeds 
in a manner similar to asymmetric stiffness effects. 

5. Double-valued responses in the form of both hard-spring 
type jump and soft-spring type jump are observed to be 
possible at low operating speeds with low auxiliary bearing 
damping or high imbalance. With large clearances or high 
bearing stiffness, the jump phenomenon may occur for 
even relatively small imbalances. 

6. The effect of friction between the shaft and the inner race 
of a rolling element auxiliary bearing on the dynamics of 
the rotor is quite small and can reasonably be neglected 
for steady-state analyses. 
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TABLE 1. FEM MODEL DATA 











(a) c = 0.0005 


(f) e - 0.0031 


0.002 

0 003 

-0.004-0.002 0.000 0.002 0.004 0.006 

X (ia) 


X (in.) 

(b) c = 0.0009 


(g) c = 0.0035 


0.002 0.003 

X (in.) 

-0.001 0.000 0.001 0.002 0.003 0.004 
X (in) 

(c) c « 0.0010 

(h) c = 0.0040 




0.002 0.003 

X (in.) 

-0.001 0.000 0.001 0.002 0.003 0.004 
X (in.) 

(d)c = 0.0015 

d 

II 

o 


300 0.002 

X (in.) 

0.004 

-0.001 0.000 0.001 0.002 0.003 0.004 
X (in.) 

(c)c = 0.0025 


(j)c = 0.0050 


0.004-0.002 0.000 0.002 0.004 0,006 
X (in.) 


-0.001 0.000 0.001 0.002 0.003 0 004 
X (in.) 


FIG. 4 IMBALANCE RESPONSES - ORBITS 
(6=0.002, A' B =0.313e+6, C b = 157.0) 










X (in.) X (in.) X (in.) X (ia) X (ia) 


(a) c = 0.0005 (0 c = 0.003 1 




(b) c = 0.0009 (g) c = 0.0035 




(c)c = 0.0010 (h)e = 0.0040 




(d)e = 0.0015 (i)e = 0.0045 




(c) c = 0.0025 (j)e = 0.0050 




FIG. 5 IMBALANCE RESPONSES - SPECTRA 
(6=0.002, A' B =0.313e+6, C B = 157.0) 
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12/2 (12 = 1400) 


12/6 (12= 1300) 


-0.004 -0 002 0.000 0,002 0.004 0.006 

-0.004 -0,002 

0.000 0.002 0.004 0.006 

X (in.) 


X (in.) 

(a) 6=0.001. e= 0.0007, C o = 157.0. 

(e) (5=0 002, 

e=0.0025, C B = 157.0. 

12/3 (12= 1400) 


12/7 (12= 1500) 


-0.004 -0.002 0.000 0.002 0.004 0.006 
X (in.) 

(b) 6=0.001, e=0.0029, C B =157.0, 


.008 -0.004 0.000 0.004 0.00 

X (in.) 

(f) 6=0.002, e=0.0032, C a = 157.0, 


0/4(12= 1700) 


0/8 ( 0 = 1200 ) 


-0.006 -0.002 0.002 0.006 
X (in.) 

-0.004 -0.002 0.000 0.002 0.004 0.006 
X (in.) 

(c) 6=0.002, e=0. 00105, C fl = 145.0, 

(g) 6=0.001, e=0.0030, C B = 157.0, 
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(d) 6=0.002, e=0. 00067, C s = 157.0, 
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(h) 6=0.002, e=0. 00395. Ca = 157.0, 


FIG. 7 SUBHARMONIC RESPONSES (ORBITS) 
(7C fl =0.313e+6) 
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FIG. 8 INFLUENCE ON CRITICAL SPEEDS 

(a) e=0.0001, C D = 157.0. /v B =0.313e+6 

(b) e=0.0001, C B = 157.0, <5=0.001 
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FIG. 9 CLEARANCE EFFECT ON JUMP PHENOMENA 
(e=0.0001, A'e=0.313e+6, Cb=300.0) 
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FIG 10 EFFECTS OF BACK-UP BEARING DAMPING 
(6=0.002, e=0.0001, A fl =0.313e+6) 




FIG. 6 CHAOTIC RESPONSES (e = 0.0023, 
6=0.002, A*B=0.313e-h6, C s = 157.0) 
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